Molecular- Spintronics: the art of driving spin through molecules 



S. SanvitcQ and A. R. Rocha 
School of Physics, Trinity College, Dublin 2, IRELAND 
(Dated: 5th February 2008) 

Spintronics is the ability of injecting, manipulating and detecting electron spins into solid state 
systems. Molecular-electronics investigates the possibility of making electronic devices using organic 
molecules. Traditionally these two burgeoning areas have lived separate lives, but recently a growing 
number of experiments have indicated a possible pathway towards their integration. This is the 
playground for molecular-spintronics, where spin-polarized currents are carried through molecules, 
and in turn they can affect the state of the molecule. We review the most recent advances in 
molecular-spintronics. In particular we discuss how a fully quantitative theory for spin-transport in 
nanostructures can offer fundamental insights into the main factors affecting spin-transport at the 
molecular level, and can help in designing novel concept devices. 
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I. INTRODUCTION 

Very few scientific discoveries have moved from an 
academic laboratory to industrial mass production as 
quickly as the giant magnetoresistance effect (GMR)i*2>, 
now exploited in any read-head for standard hard drives. 
GMR is the change of the electrical resistance of a mag- 
netic device when an external magnetic field is applied 



and it is essentially associated to a change in the mag- 
netic state of the device itself. The revolutionary sci- 
entific message revealed by the GMR effect is that the 
electron spin, as well as the electronic charge, can be 
used in electronic applications. This somehow has set a 
new paradigm. 

More recently, the electron spin has made its appear- 
ance in semiconductor physics. This new field, usually 
called spin-electronics or spintronics2i2£ has the poten- 
tial of bringing memory and logic functionalities on the 
same chip. The electron spin is the ultimate logic bit. In 
semiconductors spin preserves coherence over extremely 
long times£ and distances^, thus it offers the tantaliz- 
ing prospect of being used for quantum logic^. More- 
over all electronic ways of manipulating the spin direction 
have been proposed^. These are based on the spin-orbit 
8 interactioniS, which interestingly plays a ubiquitous role 
in semiconductor spintronics. 

One the one hand spin-orbit allows us to manipulate 
the electron spins by electric only means. It is an in- 
trinsic property of the electronic structure, and therefore 
it can be engineered by appropriate heterojunction fab- 
rication and manipulated by stress or with an external 
electric field. Importantly it can be controlled, at least 
in principle, at an extremely local level. The spin Hall 11 
effect, a solid state version of the Stern-Gerlach exper- 
iment utilizing spin-orbit interaction instead of a mag- 
netic field gradient, is a good example of all-electrical 
spin manipulation. On the other hand spin-orbit is the 
main source of spin-dephasing through the Dyakonov- 
Perel mechanism^, i.e. in semiconductors it is the main 
interaction responsible for reducing spin coherence. 

Ultimately spin looks like an attractive degree of free- 
dom to be used in logic because the energy scale relevant 
for its typical dynamics is order of magnitudes smaller 
than that involved in manipulating the electron charge 
in standard transistors. This can translate in devices 
exhibiting ultra-low power consumption and high speed. 
Moreover the sole existence of magnetic materials with 
high Curie temperature suggests the possibility of pow- 
erless non- volatility. 
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At the same time and almost in parallel there has 
been a growing interest in making electronic devices 
using organic molecules. This field, which takes the 
suggestive name of molecular-electronics 13 , aims at re- 
placing standard semiconductors with organic materials. 
These have the advantages to be manufactured with low- 
temperature low-cost chemical methods, instead of ex- 
pensive high-temperature solid-state growth (e.g. molec- 
ular beam epitaxy) and patterning (lithography) tech- 
niques. In addition the endless possibilities of chemi- 
cal synthesis and end-groups engineering give good ex- 
pectation for new concept devices. Negative differen- 
tial resistance^ and rectification^ have been already 
proved at the molecular level and prototypes of molec- 
ular transistors 17 , memories^ and logic gatesA2*A&, have 
all been demonstrated. 

It is only until recently that spin has entered the 
realm of molecular electronics. The driving idea be- 
hind the first pioneering experiment of Tsukagoshi and 
coworkers 2 ^, who injected spin polarized electrons into 
carbon nanotubes, is that spin-orbit interaction is very 
weak in carbon-based materials. This fact, in addition to 
the rather weak hyperfine interaction, suggests extremely 
long spin relaxation times, and therefore the possibil- 
ity of coherent spin propagation over large distances. A 
rather conservative estimate of the spin diffusion length 
from the Tsukagoshi's experiment indicated 130 nm as 
a lower bound of the spin-diffusion in carbon nanotubes. 
These findings have stimulated a growing activity in the 
area and several experiments dealing with molecular tun- 
neling junctions 21 , spin-transport through polymers 2 ^ 3 - 
and optical pump/probe experiments through molecular 
bridges^ have recently appeared. 

The molecular world has all the ingredients that spin- 
electronics needs. The conductivity of polymers can be 
changed by more than ten order of magnitudes 2 ^ and ele- 
mentary molecules can be designed with the desired elec- 
tronic structure. Molecules can be anchored to metals in 
numerous ways and the bonding angle can be further 
engineered by the coverage density. The spin-relaxation 
times can be extremely long and furthermore both para- 
magnetic and ferromagnetic molecules 2 ^ are available. 

Theory and modeling is a powerful engine for the de- 
velopment of this new area. At present accurate quanti- 
tative algorithms for evaluating the I-V characteristics of 
molecular devices are availablo 27 ' 28 i 29 i 30 i 31 i 32 i 33 , and they 
are revolutionarizing the world of nanoscale device simu- 
lators, as density functional theory (DFT)2i did for elec- 
tronic structure methods in the sixties. Some of these al- 
gorithms are spin polarized^LSiSL 3 ^, and therefore read- 
ily applicable to spin-transport phenomena. Certainly 
such calculations are not easy. For instance the degree 
of accuracy needed for the description of the underling 
electronic structure may go beyond what is standard in 
solid state physics 3 ^. Here in fact one needs to describe 
the metallic state of the current/ voltage electrodes, the 
molecular state of the actual device and magnetism on 
the same footing. In addition since a transport prob- 



lem is essentially a non-equilibrium problem, variational 
principles are not valid. One cannot depend on the free 
energy for atomic relaxation and the full dynamics must 
be considered 3 ^. Finally detailed information about the 
elementary excitations (phonons, spin waves etc.) and 
the exact atomic positions are essential. 

The aim of this review is to offer a complete overview 
of the fascinating field of molecular-spintronics. In par- 
ticular we will demonstrate that a quantitative theory of 
quantum transport can offer important insights and can 
be an invaluable tool for understanding complicated ex- 
periments and for novel device designing. The paper is 
organized as follows. In the first two sections we will in- 
troduce the two fields of spin- and molecular-electronics. 
Then we will introduce the main computational tools, 
and we will discuss the latest progresses with molecu- 
lar spin- valves. Finally we will overview the most recent 
and controversial findings, namely transport in molecular 
magnets, d° ferromagnetism and contact induced ferro- 
magnetism. 

II. SPIN-ELECTRONICS 
A. Transition metals and spin valves 

Magnetic transition metals and their permalloys oc- 
cupy an important place in the field of spin-electronics. 
This is essentially due to their generally high Curie tem- 
perature (for a commercially useful magnetic materials it 
must exceed ^500 °K£), and the possibility of engineer- 
ing the various magnetic properties by alloying. The fer- 
romagnetism in 3d transition metals can be understood 
by simply looking at their electronic structure. 

The nominal atomic configurations of Ni, Co and Fe 
are respectively 4s 2 3(i 8 , 4s 2 3<i 7 and 4s 2 3e? 6 . Therefore in 
forming a solid one expects the Fermi level (Ep) to be 
in a region of density of states (DOS) with dominant d 
character. Since the d shells are rather localized the DOS 
is extremely large around Ep and the material becomes 
Stoner unstable thus developing a ferromagnetic ground 
state. The band energies for the two different spin 
orientations (er =t,4) are shifted with respect to each 
other by a constant A = — , with A approximately 
1.4 eV in Fe, 1.3 eV in Co and 1.0 eV in Ni. More so- 
phisticated DFT calculations show that such picture is 
a good approximation of the real electronic structure of 
Ni, Co and Fe. 

In addition to the formation of a net magnetic moment 
a consequence of the bandstructure spin-splitting is that 
the Fermi surface for the two spin directions is rather dif- 
ferent. This difference is more pronounced in the case of 
strong fcrromagnet, where only one of the two spin-split 
d manifolds is fully occupied (majority band) while the 
other has some fractional occupation (minority band). 
An example of this situation is fee Co (the high temper- 
ature phase), whose electronic structure is presented in 
figure 2] 
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Figure 1: a) Band structure, b) density of states, and c) Fermi 
surface for fee Co. The figures al), bl) and el) refers to 
the majority spin electrons, while a2), b2) and c2) to the 
minority. The pictures a) and b) have been obtained with 
density functional theory using the code SIESTA^, and c) 
from an spd tight-binding Hamiltonian—. 



The main feature of the bandstructure of Co (and in- 
deed of the other 3c? magnetic transition metals) is the 
presence of a broad (delocalized) and only weakly spin- 
split s band, and of a narrow (localized) and largely spin- 
split d band. The former has almost free electron-like 
character for energies both below (note the parabolic be- 
havior of the bands around V for E ~ -9 eV) and above 
E-p. In contrast the latter is only about 5 eV wide and 
cuts close to the Fermi level. Because of the spin-splitting 
of the d manifold, A, the majority spin band has an al- 
most spherical Fermi surface (see figure ^1), while the 
minority one has a rather complicated structure, mostly 
arising from the d manifold (see figure |TJ;2). The differ- 
ent structure of the Fermi surface for majority and mi- 
nority spins and the fact that these differences arise from 
a different orbital character are the main ingredients for 
understanding the transport properties of magnetic tran- 
sition metal heterostructures. 



Importantly the effects of alloying can be understood 
in the context of a rigid-band model, by shifting the 
Fermi level according to the valence of the dopanl^Mi. 
Finally it is worth mentioning that there are materials 
that at the Fermi level present a finite DOS for one spin 
specie and a gap for the other. These are known as half- 
metals^ and are probably among the best candidates as 
materials for future magneto-electronics devices. 



B. Spin valves 

The prototype of all spin devices is the spin-valve. 
This is formed by two magnetic layers (normally tran- 
sition metals) separated by a non-magnetic spacer (ei- 
ther metal or insulator). One of the two magnetic layers 
is free to rotate in tiny magnetic fields, while the other 
usually is pinned by exchange coupling with a antiferro- 
magnet or by strong magnetic anisotropy. The current 
passing through a spin-valve depends over the mutual 
orientation of the two magnetic layers and it is typically 
higher for a parallel alignment (PA) than for an antipar- 
allel (AA) . Thus a spin- valve behaves essentially as a spin 
polarizer/analyzer device. The quantity that defines the 
effectiveness of the spin-filtering effect is the GMR ratio 
]?gmr defined as ("optimistic definition") 



-Rgmr 



h 



PA 



(1) 



An alternative definition ( "pessimistic definition" ) using 
Jpa + Iaa as normalization is sometime used. 

An intuitive understanding of the spin-filtering pro- 
duced by a spin-valve can be obtained by looking at the 
Fermi level lineup of the materials forming the device. 
Let us consider for example a Ni/Cu/Ni spin- valve (see 
figure^ , and assume that the two spin-bands do not mix. 
This is the two spin fluid approximation, which is valid 
in the case of weak spin-orbit scattering and collinear 
magnetism^. 

In the AA a spin electron propagates in the majority 
spin band in one layer and in the minority band in the 
other. Consequently electrons always travel across the 
Fermi surfaces of Cu and of both the spin-bands of Ni. 
In contrast in the PA the two spin currents are rather 
different. The majority current is made from electrons 
that have traveled within the Fermi surfaces of Cu and 
that of the majority spin of Ni, while the down spin cur- 
rent from electrons that have traveled within the Fermi 
surfaces of Cu and that of the minority spin of Ni. This 
leads to two different current paths for the AA and the 
PA. Since the two spin currents add to form the total cur- 
rent and since generally the resistances of majority and 
minority electrons are different in a magnetic transition 
metal, the current passing through the PA and AA con- 
figurations are different. Importantly the larger is the 
mismatch between the two spin currents, the larger in 
the GMR ratio. 
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Figure 2: Magnetoresistance mechanism in a Ni/Cu/Ni spin 
valve in the two spin fluid approximation. In the parallel case 
a) the resistance of the majority spin channel is low since there 
is good match of the majority Fermi surfaces across the entire 
device. In contrast in the antiparallel state b) the alignment of 
the Fermi surface is such to have one high resistance interface 
for both the spin channels. This is given by the interface 
between the minority spin band of Ni and that of Cu. The 
formation of high resistance channels for both spins in the AA 
is responsible for the GMR effect. 



C. Spin polarization of a device 

An important question is how to quantify the relative 
difference between the two spin currents in a magnetic 
device and how to relate this properties to the elemen- 
tary electronic structure of the materials forming the de- 
vice. We then define the spin-polarization P of a mate- 
rial/device as 



P 



h 



J T + /,. ' 



(2) 



where I a is the spin-er contribution to the current. l a and 
P are not directly observable and must be calculated or 
inferred from indirect measurements. Unfortunately the 
way to relate the spin-current I a to the electronic struc- 
ture of a material is not uniquely defined and depends on 
the particular experiment carried out. 

As brilliantly pointed out by Mazini4, the relation be- 



tween the spin-polarization of a magnetic material and 
its electronic structure depends critically on the trans- 
port regime that one is considering (ballistic, diffusive, 
tunneling ..). As a first approximation the current I is 
simply proportional to Npv F , where N-p and vp are the 
DOS at the Fermi level and the Fermi velocity respec- 
tively. Different transport regimes weight the contribu- 
tion of the Fermi velocity differently, and one has n = 2 
for diffusive transport, n = 1 for ballistic transport and 
n = for tunneling. 

Therefore the spin-polarization P becomes 



N F (v F ) n - N F (v F y 



(3) 



Typical values of P n for several magnetic metals are re- 
ported in table U 



Pn (%) n=2 n=l n=0 



Fe 


20 


30 


60 


Ni 





-49 


-82 


Cr0 2 


100 


100 


100 


Lao.67Ca .33Mn0 3 


92 


76 


36 


Tl 2 Mn 2 7 


-71 


-5 


66 



Table I: Spin-polarization of typical magnetic metals accord- 
ing to the various definitions given in the text. The data 
are taken from literature as follows: Ni and Fa—, CrOa—, 
Lao.erCao.ssMnOs 46 and Tl 2 Mn 2 OA 

Importantly the spin-polarization of a device can be 
different from that of the materials forming it. This 
is connected to the fact that the bonding at the inter- 
face between two different materials can be strongly spin- 
selective. For instance if the bonding between two ma- 
terials has mainly s-character, then one expects strong 
scattering for d-like electrons. As a consequence the 
spin-polarization of the current will be determined by 
the spin-polarization of the almost free s-electrons. This 
is usually much smaller than that of the (i-electrons and 
it may even have the opposite sign. For instance it is 
demonstrated that the sign of the magnetoresistance of 
a magnetic tunneling junction can be altered by simply 
replacing the insulator forming the barrier—. 



III. MOLECULAR-ELECTRONICS 

A. Electrons transport through molecules 

Transport through molecules and in general through 
low dimensional objects is somehow different than that 
in standard metals or semiconductors. This is substan- 
tially due to the collapse of the Fermi surface into a sin- 
gle energy level (the highest occupied molecular orbital 
- HOMO). The nature and lineup of the HOMO with 



5 



the Fermi energy of the current /voltage probes deter- 
mine most of the transport properties. Let us consider 
the simple case of a two probe device. Following a sim- 
ple model proposed by Datta^ the typical energy level 
lineup is schematically presented in figure [3J 



a) 



b) 
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Figure 3: Energy level lineup between a molecule and two 
current voltage probes. In the weak coupling limit a) the 
molecule is characterized by the ionization potential (IP) and 
the electron affinity (EA) , which line up with the metal Fermi 
energy (WF is the work function of the contact). In the case 
of strong coupling b) between the molecule and the leads the 
molecular levels shift and broaden. It is then more appro- 
priate to discuss transport in terms of HOMO and LUMO 
states. 

In absence of any coupling (figure Ot) both the energy 
levels of the molecule and the Fermi level of the elec- 
trodes will align with a common vacuum level. In this 
case the system is characterized by the work function of 
the electrodes and both the ionization potential (IP) and 
the electron affinity (EA) of the molecule. In this setup 
the molecule can exchange electrons with the electrodes 
only if the typical temperature is comparable to cither 
IP-WF or WF-EA, a condition which is normally not 
satisfied. This guarantees local charge neutrality of the 
whole system and integer occupation of the molecule. 

In contrast, the interaction between the molecular lev- 
els and the extended wave-functions of the metallic con- 
tacts has the effect of broadening and shifting the molecu- 
lar levels. In the extreme limit of large coupling extended 
states spanning through the entire system (electrode plus 
molecule) can develop and the molecular device will be- 
have as a good metal. In this limit the molecular lev- 
els cannot be associated any longer to the elementary 
removal energies of the isolated molecule and a descrip- 
tion in terms of fractionally occupied HOMO and LUMO 
(lowest unoccupied molecular orbital) is more appropri- 
ate. 

The transition from integer to fractional occupation of 
the molecule somehow depends on the typical molecular 
charging energy (say the EA) compared to the hopping 
integral T between the molecule and the contacts (F/K 
is the escape rate from the molecule to the contacts). 
One has integer occupation if EA ^> T and metallic-like 
behaviour when T ^> EA. This is essentially the same 
physics leading to Mott metal-insulator transition in solid 
state. An electronic structure theory capable of explor- 
ing on the same footing all the intermediate situations 
between the strong and weak coupling limit is still not 



available unless at prohibitive computational costs^ 5 .. 

The effect of an applied bias V is that of shifting the 
chemical potentials of the two current/ voltage probes rel- 
ative to each other by eV, with e the electronic charge. 
As a rule of thumbs current will flow whenever a molecu- 
lar level (either the HOMO or the LUMO) is positioned 
within such a bias window. The appearing of molecu- 
lar levels in the bias window when the potential is in- 
creased typically leads to changes in the slope of the I-V 
characteristics, in steps in the differential conductance 
dI/dV(V) and in peaks in its derivative d 2 I/dV 2 (V). 
This means that fingerprints of the molecular spectrum 
can be found in the measurement of its electrical proper- 
ties. The same is true for the molecular elementary ex- 
citations, and peaks in d 2 I/dV 2 (V) can be found in cor- 
respondence of the energy of relevant phonon modest. 



B. The bonding with the contacts 

One of the fundamental aspects of molecular electron- 
ics is that the bonding between a molecule and the cur- 
rent/voltage probes can be engineered to a degree usually 
superior to that achievable in conventional inorganic het- 
erostructures. This can dramatically change the current 
flowing through a device. Consider for instance the sim- 
ple case of an atomic gold chain, described by s orbitals 
only, sandwiching a 7r-bonded molecule (see figure 0J. 




P. P. 




Figure 4: Au atomic chain sandwiching a 7r-bonded molecule 
(say an S2 molecule), a) the molecule is aligned with the Au 
chain and the transmission is suppressed because the matrix 
element (s\H\p x ) vanishes. In contrast when the molecule 
forms some angle 6 with the Au chain b) then a component 
of the hopping integral along the bond develops and current 
can flow. 

Let us assume that the relevant molecular state (the 
one close to the Fermi level of the gold chain) is formed 
by p x orbitals, i.e. those perpendicular to the chain axis. 
When the molecule is positioned along the axis of the 
chain the hopping integral between the molecule and the 
chain (s\H\p x ) vanishes regardless of the separation be- 
tween the two. This is simply the result of the particu- 
lar symmetry of the problem since s and p x orbitals do 
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not share the same angular momentum about the bond 
axis 51 . As a consequence the current is identically zero. 

In contrast if the molecule is not coaxial to the chain 
(figure EJ)), then there is a component of the p x orbital 
along the bond axis and the hopping integral becomes 
"fspa sin (9, where "f S pa is the spa hopping integral and 
9 the bond angle. This dependence of the bonding on 
the bond orientation may have dramatic consequences 
on the I-V characteristics. In figure we present the 
I-V curve for the system of figure 01 calculated with a 
self-consistent tight-binding model where it is assumed a 
linear dependence of the on-site energies over the orbital 
occupation (see reference [3j| for details). 
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Figure 5: I-V characteristics for a S2 molecule sandwiched 
between two semi-infinite gold chains (see figure The 
curves are calculated using the tight-binding based non- 
equilibrium Green function method of reference— with the 
following parameters: £a u =-5.9 eV, es=-6.15eV eV, j S sa=- 
3.0 eV, 7 spCT =1.52 eV, 7 Ppw =-0.98 eV, t/ Au =-6.7 eV and U s =- 
6.15 eV. e is the on-site energy, 7 hopping integral and U the 
charging energy. 

From figure 03 it is clear that the I-V characteristics 
of a molecule can be largely engineered by simply chang- 
ing the details of the bonding with the electrodes. This 
is considerably more complicated in extended interfaces 
(for instance between two metals), since disorder, inter- 
diffusion and roughness have the effect of averaging out 
the atomistic details of the bonding. It is important to 
remark that even when molecular layers are grown a good 
level of tuning of the bonding properties still exists. For 
instance the bonding site and the bond angle usually de- 
pend on the layer density (coverage) 5 ^, and these can be 
further tuned by changing the end groups. 



C. Why spins and molecules? 

What are the advantages of using molecules instead of 
inorganic materials for performing spin-physics? These 
are essentially two. On the one hand there are intrin- 
sic molecular properties and in particular the weak spin- 
orbit and hyperfine interactions. On the other hand there 



are the properties connected to the formation of inter- 
faces between magnetic metals and molecules. We will 
return to the interfacial properties in the next sections, 
here we focus our attention only on the intrinsic aspects. 

Spin-orbit interaction is a relativistic effect which cou- 
ples the electron spin S with its angular momentum L. 
The spin-orbit Hamiltonian in general can be written as 
^so = ^SoS" ' L, where Vso(f) is a term which contains 
the gradient of the electrostatic potential. Although it 
is rather intuitive to realize that the strength of this in- 
teraction grows with the atomic number Z (it is propor- 
tional to Z 4 ), its actual value in the solid state depends 
on various factors such as the crystal symmetry and the 
material composition. Importantly the spin-orbit effect 
is responsible for spin-precession and the loss of spin- 
coherence. In organic materials usually the spin-orbit in- 
teraction is rather small. This is mostly due to the small 
atomic number of carbon. In table ^ we compare the 
spin-orbit splitting Ago of the valence band of several 
semiconductors 5 - with that of carbon diamond 54 . The 



Aso (meV) 


Si 


44 


Ge 


290 


GaAs 


340 


AlAs 


280 


InAs 


380 


GaP 


80 


InP 


111 


GaSb 


750 


AlSb 


670 


InSb 


980 


C 


13 



Table II: Valence band spin-orbit splitting for various semi- 
conductors. 

table shows that in carbon the spin-split of the valence 
band is approximately one order of magnitude smaller 
than in ordinary III-V or group IV semiconductors and 
one should expect a considerably longer spin- lifetime 6-7 . 

Another important interaction, which generally leads 
to spin-decoherence, is the hyperfine interaction between 
electron and nuclear spins. This has the form 

-ffhyp = ^hyp S • Sn , (4) 

where s and Sn are respectively the electronic and nu- 
clear spin, and Ah yp is the hyperfine coupling strength. 
Similarly to the case of spin-orbit, also hyperfine interac- 
tion is a source of spin de-coherence^ 5 -, since the random 
flipping of a nuclear spin can cause that of an electron 
spin. However, in III-V semiconductors hyperfine inter- 
action was also proved to be a tool for controlling nuclear 
spins via optically polarized electron spins£&. In organic 
materials usually the hyperfine interaction is weak. The 
main reason for this is that most of the molecules used 
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for spin-transport are 7r-conjugate molecules where the 
transport is mostly through molecular states localized 
over the carbon atoms. Carbon, in its most abundant 
isotopic form, 12 C, has nuclear spin 5n=0, and therefore 
is not hyperfine active. Moreover the 7r-states are usually 
delocalized and -ffhyp can be anyway rather small. In ta- 
ble lllll we report the value of the nuclear spin for various 
atomic species with their relative isotopic abundance. 



Isotope IA (%) 


S'n 


X H 


99.98 


1/2 


2 H 


0.02 


1 


12 C 


98.93 





13 C 


1.1 


1/2 


14 N 


99.632 


2 


15 N 


0.368 


1/2 


16 o 


99.757 





18 o 


0.205 





19 F 


100 


1/2 


69 Ga 


60.108 


3/2 


71 Ga 


39.892 


3/2 


75 As 


100 


3/2 


28 Si 


92.2297 





29 Si 


4.6832 


1/2 


30 Si 


3.0872 






Table III: Nuclear spin for elements present in typical organic 
molecules and in both Si and GaAs. Here we report the nu- 
clear spin for the most abundant isotopes, together with their 
relative isotopic abundance (IA). 

Estimates of the spin- lifetime of organic materials from 
transport experiments are at the moment only a few. 
Moreover in most cases these are extracted from spin- 
valves measurements by fitting to the Jullier's formula 57 . 
This procedures does not distinguish the source of spin- 
flip, which may not be located inside the molecule, but 
at the interface with the magnetic electrodes. Therefore 
these measurements are likely to offer a conservative esti- 
mate of the spin-lifetime. Nevertheless the values of the 
spin diffusion length reported in the literature are rather 
encouraging for carbon nanotubes (130nm)SS, polymers 
(200nm)2i, or Alq 3 molecules (5nm)£&. 

Finally a conclusive note needs to be made on con- 
ducting polymers. These are extremely attractive ma- 
terials since their electrical conductivity can be changed 
by over twelve orders of magnitude with the chemical 
or electrochemical introduction of various counterions 59 . 
Spin-dynamics in polymers has been extensively stud- 
ies with EPS spectroscopy^ and it is largely dominated 
by the presence of paramagnetic centers in the form of 
free radicals, ion-radicals, molecules in triplet states and 
transition metal complexes. Even more interesting is the 
fact that the elementary excitations leading to electron 
transport are not band-like but usually involve lattice 
vibrations, and most importantly some of them are spin- 



polarized. 

For instance in the irans-isomer of polycetylene the Su- 
Schrieffer-Heeger theory^ predicts a soliton-like trans- 
port mechanism. This has a peculiar spin-charge rela- 
tionship, since a neutral soliton corresponds to a radi- 
cal with spin 1/2, while both negatively and positively 
charge solitons are spinless and diamagnetic. The study 
of spin-transport in devices made by transition metals in 
contact with such polymers is potentially extremely in- 
teresting. This is because of the proximity of two funda- 
mentally different ground states, one electronically corre- 
lated (the magnetic material) , and one with strong corre- 
lation between electronic and vibrational degrees of free- 
dom (the polymer). The investigation, both experimen- 
tally and theoretically, of these combined systems is in 
its infancy 61 . 



IV. QUANTITATIVE TRANSPORT THEORY 
A. Simple Model 

Modern theory of quantum transport is based on scat- 
tering theory in conjunction with accurate electronic 
structure methods. Although this approach has recently 
come to question, in particular in the case of electrons 
interacting beyond the mean-field level 6 ^, it still remains 
the most versatile and scalable available. In addition its 
foundations are extremely intuitive and simple. Let us 
start our discussion by presenting a simple model, first 
introduced by Datta 49 , which already contains all the el- 
ements of a more formal and accurate theory. 

Consider a given molecule attach to two cur- 
rent/voltage probes kept at two different chemical po- 
tentials fi a , with a=L (left), R (right). The leads 
are assumed featureless, which is with a constant DOS. 
The molecule is described by an energy level e (say the 
HOMO), coupled to the current/ voltage probes by the 
hopping integrals t a (a=L, R). The density of states as- 
sociated to such a state is 



D e (E) 



t/2-K 



{E-ef + {t/2f 



(5) 



where the level lifetime Ti/tvs determined by the coupling 
with the leads only t = £l+£r- Both the current / flowing 
through the molecular state and state occupation N can 
be determined by balancing the in-going and out-going 
fluxes and read 

/ + OG ii 
dED £ (E)±f[f L (E)-f R (E)}, (6) 



and 



N 



(7) 



/ di?W L/L(ig)+WR(ig) 

J — oo * 

where f a {E) = 1+e (,- l a )/k B T 1S the Fermi distribution of 
the contact a. From the equations above it is clear that 
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current will flow only (at least for small broadening and 
low temperature) if the energy level is in between the 
chemical potentials of the two leads. With no external 
bias these are identical, but when a potential V is applied 
then /il — Pn = eV and current will flow. 

In general e depends on the details of the electronic 
structure of the system. As a simple approximation we 
may assume it is a function of the level occupation only 
e = e(N). This suggests a simple self-consistent proce- 
dure where e(N) and the equation JJJ) are solved itera- 
tively before the current is evaluated with the (JBJ. 

Spin can be easily introduced in this simple model in 
the two spin fluid approximation. Since the two spin 
channels do not mix the equations © and JJJ) can be 
replaced by two pairs of equations for the spin-resolved 
molecular level occupation N a and the spin-current I a . 
In a similar way the molecular level can also be spin- 
polarized e — > e CT and it may depend of the spin den- 
sity instead of the density only e a (N^ , N^-). Importantly 
in the case where the current/ voltage probes are ferro- 
magnetic spin-degeneracy is lifted by introducing spin- 
dependent hopping integrals t" a between the leads and 
the molecule. These however capture the fact that ma- 
jority and minority electrons couple to the molecule in 
a different way, but not that the DOS for the different 
spin directions in a ferromagnet is different. For this 
last feature a more detailed description of the electronic 
structure of the leads is needed. 



B. Ab initio methods 

The non-equilibrium Green's function (NEGF) method 
is by far the most used among all the quantum transport 
schemes. Although it is based on very rigorous ground 63 
it can be understood as the natural extension of the toy- 
model discussed in the previous section. The general idea 
is to divide a two probe device into three regions: two 
current/ voltage probes and a scattering region. The cri- 
terion for this fragmentation is that the scattering region 
is the portion of the device where the potential drops. Its 
boundaries are defined by the condition that the charge 
density matches exactly that of the bulk material forming 
the leads (see reference^ 7 - for a more detailed description) . 

Let us assume that the problem can be formulated over 
some sort of localized basis set, and therefore the Hamil- 
tonian for the whole system (scattering region plus leads) 
is simply an infinite hermitian matrix. The central quan- 
tity, which replaces the simple density of states D t (E), is 
the non-equilibrium Green's function for the scattering 
region G(E) 

G(E) = limlOE + ^-ffs-EL-SR,]- 1 . (8) 
77— >o 

This is the Green's function associated to the Hamilto- 
nian TJs+El+Er, which is composed by the Hamiltonian 
of the scattering region Hs and the self-energies El and 
Sr. 



Note that if one assumes that Hs, El and Er are just 
C-numbers, then i[G(E) - G*{E)] is D e (E) for e -> H s 
(real) and t a /2 — > E Q (purely imaginary). Thus the self- 
energies can be associated with the interaction between 
the scattering region and the current /voltage probes. In 
practice they can be written as El = H^g^H^s and 
Er = iJRsgR.if^g, with H a s the coupling matrix between 
the leads a and the scattering region. g a are the surface 
Green's functions for the leads, i.e. the Green's func- 
tion for a semi-infinite lead evaluated at the termination 
plane^i. Thus the self-energies contain information on 
both the coupling between the scattering region and the 
electrodes and on the electronic structure of the leads 
themselves. 

The crucial point is that both the two-probe current / 
and the density matrix associated to the scattering region 
p can be obtained from the Green's function G(E) 

I=^J dE Tr{GT L G^T R }[f L (E) - ME)} , (9) 
and 

p= hj dE g[iy/l + Tr ^ e ^ > ( 10 ) 

where G = G(E), T = T(E) we have now introduced the 
broadening matrices 

r a = *[s a -st]. (11) 

Equations (JSJl and 10 allow us to calculate the cur- 
rent once the Hamiltonian for the scattering region and 
the self-energies are given. Unfortunately these are not 
known a priori, since they require the evaluation of the 
electronic structure of an infinite non-periodic system. 
However their calculation do not require the actual solu- 
tion of this open system and the notion of locality (the 
same that allows us to separate the leads from the scat- 
tering region) can be efficiently used. The self-energies in 
fact involve only bulk quantities, i.e. quantities which can 
be evaluated from a bulk calculation for a periodic sys- 
tem. In addition a self-consistent scheme can be designed 
for evaluating H§. The crucial point here is to assume 
that the whole electronic structure can be described by 
a single-particle theory, in such a way that Hs depends 
only on the density matrix Hs = Hs[p\- This is for in- 
stance the case of DFT or Hartree-Foch methods. Hence 
H s [p] and the equations © and ljTU|l can be iterated self- 
consistently until convergence, and then the current can 
be evaluated with equation This scheme, with differ- 
ences concerning the specific numerical implementation, 
is used by most of the DFT-based NEGF packages avail- 
able at present 27 . 28 i 29 i 30 i 31 . 32 i 33 . 

Also in this case the addition of spin polarization does 
not bring any fundamental changes in the formalism. In 
the simple case of collinear-spins, then all the quantities 
introduced (Green's functions, charge density, Hamilto- 
nian, self-energies ...) becomes block diagonal matrices 
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in spin space. The Hamiltonian now depends on both 
the charge density and the magnetization as in standard 
spin-polarized mean field electronic structure methods. 
Importantly the current is still carried in parallel by the 
two spin species. In the more complicated case of non- 
collinear spin (or if spin-orbit is present) then the off- 
diagonal blocks to not vanish and the current cannot be 
broken down into the majority and minority contribu- 
tions. 

Note that computationally the introduction of mag- 
netism usually makes the calculation more complicated. 
First one needs to consider matrices larger than the 
unpolarized case. Secondly and most importantly the 
Hamiltonian matrix becomes considerably more sparse. 
In a ferromagnet in fact, for instance in a transition 
metal, the delocalized s-electrons responsible for most 
of the electron transport coexist with the tightly bound 
d-electrons, which provide the local magnetic moment. 
This means that most of the matrix elements between 
d orbitals located on atoms far from each other vanish, 
making the Hamiltonian more sparse than that of stan- 
dard free-electron like metals (for instance Au). In the 
case of sparse matrices standard recursive methods for 
evaluating the self-energies become numerically unstable 
if not prohibitive, and more sophisticated schemes are 
neededSLSa. 



O O O 

O O O 

O O 

O O O 

o o o 




o o o 

QOO 
3 O O 

o o o 

o o o 



Figure 6: Setup for early molecular spin- valve calculations. 
A molecule (black circles) is sandwiched between two mag- 
netic ions (blue circles), which are then contacted by two non- 
magnetic current /voltage electrodes. The red arrows indicate 
the local magnetic moments. 



An interesting alternative, that still avoids the con- 
struction of magnetic leads, was proposed by Wei and 
co-workers^, who used Al leads locally immersed in a 
magnetic field. In this case the spin-polarization of the 
leads is achieved through simple Zeeman splitting and the 
direction of the local magnetic fields at the two contacts 
replaces that of the magnetic moments of a ferromag- 
net. In this case the DOS of the leads is spin-polarized, 
however the orbital nature of the bonding is identical for 
both majority and minority spins. Importantly most of 
these early calculations were ahead of experiments and 
largely inspired them. 



MOLECULAR SPIN- VALVES 



As pointed out in the introduction spin-valves are the 
prototypical spin-devices and therefore are the most stud- 
ied architectures for spin-transport through molecules. 
The first prediction of a GMR-like effect in molecules 
is from Emberly and Kirczenow, who investigated the 
transport through a 1,4-benzenedithiolate molecule at- 
tached to Ni electrodes. Because of the problems of deal- 
ing with magnetic current /voltage probes mentioned in 
the previous section, this work and most of the early 
calculations^ were limited to empirical models for the 
electronic structure. 

A first step in the direction of using ab initio elec- 
tronic structure methods was suggested by Pati and co- 
workersSSiSLS, who considered a setup in which a given 
molecule is sandwiched between two magnetic ions, or 
clusters, which then are attached to non-magnetic elec- 
trodes (see figure El- 

This setup removes the problems connected with con- 
structing the self-energies for magnetic materials and a 
GMR-like effect is produced by the magnetic ions, which 
in turn are exchanged coupled through the molecule. 
Clearly the scheme is highly idealized. It essentially 
describes the transport through a magnetic molecule 
(molecule plus magnetic ions) from non-magnetic leads, 
and not that of a magnetic spin-valve. Importantly it 
does not account for both the spin-polarized DOS of the 
contacts and the accurate orbital character of the bond- 
ing between the molecule and the magnetic surfaces. 



A. Metallic and tunneling junctions 

The study of molecular transport with ferromagnetic 
contacts from first principles has a young history since 
only recently algorithms stable enough to deal with ex- 
tremely sparse matrices and thus with magnetic leads 
were made availabl o 27 i 28 . Here we describe the results 
obtained with the code Smeago? 7 - 28 for both insulating 
and tunneling molecules. 

As an example of different transport regimes, we 
consider spin-valves made from Ni leads and a 
molecular spacer which is either [8]-alkane-dithiolate 
(octane-dithiolate) or I,4-[3]-phenyl-dithiolate (tricene- 
dithiolate). A schematic DOS and the charge density 
isosurfaces of the HOMO and LUMO states for the iso- 
lated molecules are presented in figures and El 

The two molecules present rather different characteris- 
tics. The HOMO-LUMO gapZ&is about 2.5 eV for tricene 
and almost double (5 eV) for octance. In addition while 
in tricene the charge density of the first two HOMO lev- 
els and the LUMO is extended over the whole molecule, 
in octane this is predominantly localized around the S 
atoms of the thiol groups. We then expect that the 
octane and the tricene will form respectively TMR and 
GMR devices, as actually found in our calculations 2 ^. 

Let us consider octane first. The zero-bias transmis- 
sion coefficient of the Ni/octane/Ni junction presents 
a sharp peak at E-p that scales exponentially with the 



10 




Figure 7: [8]-alkane (octane) molecule: DOS and charge den- 
sity isosurface plots for the relevant molecular states of the 
isolated molecule. Ef denotes the position of the Fermi level 
for the isolated molecule. 
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Figure 9: ln[T(_Ep)] as a function of the number of C atom 
for [n]-alkane-dithiolate attached to Ni leads. The black cir- 
cles correspond to our calculated values, and the solid line 
is our best linear fit. Here the spin- valve is in the parallel 
configuration. 
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Figure 8: 1,4- [3] -phenyl (tricene) molecule: DOS and charge 
density isosurface plots for the relevant molecular states of 
the isolated molecule. Ef denotes the position of the Fermi 
level for the isolated molecule. 



number of alkane groups T oc e~^ n ((3 ~ 0.88). This 
is demonstrated in figure [5] for the parallel configura- 
tion and confirms that device is in a tunneling regime. 
Note that the method is accurate down to a conduc- 
tance of ~ 50pS. Interestingly the exponent is similar 
to that found for the same molecule attached to gold 
(111) surfaces^. The coupling between the thiol groups 
and the electrodes is strong and it gives rise to a small 
spin-polarization of the two S atoms. However since both 
the HOMO and LUMO states are strongly localized at 
the thiol groups, at the Fermi level there is no molecular 
state extending through the entire structure. Hence the 
Ni/octane/Ni junction presents the features of a TMR 
spin- valve. 

The case of l,4-[n]-phenyl-dithiolate is different, in par- 
ticular we find that the current does not scale sensibly 
with the number of phenyl groups in the molecule. In fig- 



ure El w e present the transmission coefficient as a func- 
tion of energy of l,4-[n]-phenyl-dithiolate for 1, 3 and 4 
phenyl rings when the magnetization vectors of the leads 
are parallel to each other. It is clear that, although there 
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Figure 10: Transmission coefficient as a function of energy 
for l,4-[n]-phenyl-dithiolate molecules with n=l, 3, and 4. 
Here the spin- valve is in the parallel configuration. A similar 
scaling of the transmission coefficient as a function of n is 
found for the antiparallel configuration. 

is a reduction of the transmission coefficient as a func- 
tion of the number of rings, this remains close to unity for 
most of the energy range, and certainly there is not an ex- 
ponential decay. Therefore in this class of molecules the 
transport seems to be appropriately described by a co- 
herent resonant tunneling mechanism through extended 
molecular states. 

This picture is enforced by the fact that the zero- 
bias transmission coefficient approaches unity for ener- 
gies close to the leads Fermi level. In addition, from 
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the study of the evolution of the orbital resolved density 
of states as a function of the distance between the thiol 
group and the electrodes^ 7 - we identify such a resonant 
state as the HOMO state of tricene. However it is worth 
mentioning that this appears rather broad and spin-split, 
because of the strong coupling with the d orbitals of the 
leads. In conclusion all this suggests that a Ni/tricene/Ni 
spin-valve behaves as metallic spin-valve. 

The I-V characteristics of both the molecules are 
strongly non-linear with the bias and consequently also 
the GMR ratio suffers this non-linearity. In figures El 
andll2lwe present the I-V curves, the GMR ratio and the 
zero-bias transmission coefficient for the two molecules. 
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Figure 11: a) I-V characteristic, and zero bias transmission 
coefficients for the b) parallel and c) antiparallel configura- 
tion of an octane-based Ni spin-valve. In the antiparallel case 
the transmission coefficient is identical for both the spin di- 
rections. In the inset we present the corresponding MR ratio. 
_Ef is the position of the Fermi level of the Ni leads. From 
reference-^. 

Again consider octane first. Here the transmission co- 
efficient at zero bias is dominated by a number of sharp 
peaks in the parallel configuration, which get consider- 
ably suppressed in the antiparallel one. In particular a 
minority peak appears at the Fermi level and it is mostly 
responsible for the low-bias conductivity. Such sharp 
peaks in the transmission coefficient are usually a signa- 
ture of resonant states at the interface between the leads 
and the insulating media. At resonance they can carry 
a considerable current, however small bias and disorder 
are usually rather effective in suppressing their contribu- 
tion to the current. Evidence for these surface states has 
been already provided for conventional magnetic tunnel- 
ing junctions both experimentally 7 ^ and theoretically 73 . 

The I-V characteristic is rather linear for the antipar- 
allel configuration but presents a non-trivial slope for the 
parallel case. This gives rise to a non-monotonic depen- 
dance of the GMR ratio over the bias. Importantly both 
the layer resistance and the TMR ratio are in the same 
range as in recent experiments on octane-based Ni spin- 
valvesSi. However, a direct comparison with experiments 



40.0 

^20.0 
0.0 
"-20.0 
-40.0 

1.5 

@ 1 
H 0.5 



g 1 
H0.5 





a) 



P configuration 
AP Configuration 



-l 





V (Volt) 



F ' } 

lA i 


i 




: c) 


1 

,1k Ma 


— — Majority 


Minority ^^~\^~^ 





(E-E ) [ eV ] 



Figure 12: a) I-V characteristic, and zero bias transmission 
coefficients for the b) parallel and c) antiparallel configura- 
tion of a tricene-based Ni spin-valve. In the antiparallel case 
the transmission coefficient is identical for both the spin di- 
rections. In the inset we present the corresponding MR ratio. 
Ef is the position of the Fermi level of the Ni leads. From 
reference—. 



is difficult since the actual number of molecules bridging 
the two electrodes is not known with precision. Moreover 
a degradation of the GMR signal due to spin-flip and 
electron-phonon scattering, misalignment of the magne- 
tization of the contacts and current shortcut through 
highly conductive pin-holes, can drastically reduce Rqmk 
in actual samples. 

In contrast in the case of metallic (tricene) junctions 
the transmission coefficient at zero bias presents a much 
higher transmission and is a rather smooth function of 
the energy. In the parallel case most of the transmission 
at Ep is due to the majority spin, while the contribution 
of the minority is significant only for energies at about 
200 mcV above Ep. In the antiparallel configuration such 
transmission at the Fermi level is strongly suppressed and 
the resulting current is considerably lower. This produces 
an extremely large GMR ratio for small bias, exceeding 
600%. 

An interesting feature is that, in first approximation, 
the transmission coefficient at zero bias for the antipar- 
allel state appears to be a convolution of those for the 
majority and minority spin in the parallel case. This find- 
ing can be qualitatively understood in terms of transport 
through a single molecular state (see figure lT5|l . Let (E) 
be the majority spin hopping integral from one of the 
leads to the molecular state, and t^(E) the same quan- 
tity for the minority spins. Then, neglecting multiple 
scattering (i.e. all interference effects), the total trans- 
mission coefficients of the entire spin- valve in the parallel 
state can be written (E) = (t r ) 2 and T U (E) = (t 1 ) 2 
respectively for the majority and minority spins. Simi- 
larly the transmission in the anti-parallel configuration 
is T^(E) = T^(E) = tUK Thus T^(E) is a convolu- 
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tion of the transmission coefficients for the parallel case 





Figure 13: Scheme of the spin-transport mechanism through 
a single molecular state, (E) (t (E)) is the majority (mi- 
nority) spin hopping integral from one of the leads to the 
molecular state. Neglecting quantum interference, in the 
parallel case (a) the total transmission coefficient is simply 
T = {t r f + (t 1 ) 2 , while in the antiparallel (b) T = 2th l . 
Note that if either or tr vanishes, the current in the an- 
tiparallel configuration will also vanish (infinite GMR). 

An extreme case is when only one spin couples to the 
molecular state. Then the total transmission in the an- 
tiparallel case is identically zero since either t* or van- 
ishes. This is the most desirable situation in real devices 
since, in principle, an infinite -Rgmr can be obtained. 
Note that in this situation the system leads+molecule 
behaves as a half-metal although the two materials form- 
ing the device are not half-metals themselves. An even 
more extreme situation is when for a particular energy 
window the transport is through two distinct molecular 
states, which are respectively coupled to the majority and 
minority spin only. This may happen for instance due a 
particular symmetry of the molecular anchoring groups. 
Then in this energy window one will find T^(E) ^ 0, 
T ll (E) ^ but T n {E) = (see figureEJ). 

The fact that in good approximation <x y/T^T^ 
can be used for enhancing the GMR ratio. Consider 
for instance the case when TJ T 3> T*+. Then we have 
T P « Tfl and T A p ~ VTTTTii for the total transmission 
coefficients respectively of the parallel and antiparallel 
configurations. Clearly a reduction of will produce 
a considerable reduction in the transmission of the an- 
tiparallel alignment, leaving almost unchanged that of 
the parallel one. We have explored this avenueSI by re- 
placing the thiol group with either a Se or a Te atoms. 
These provide a rather strong bond, although generally 
the bond length is increased due to the larger atomic ra- 
dius of the anchoring atom. The transmission, I-V char- 
acteristics and GMR for 1,4-benzene anchored to Ni via 
S, Se or Te are presented in figures IT51 ITrjl and IT71 Clearly 
the pictures show a rather dramatic increase of the GMR 
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Figure 14: Scheme of the spin-transport mechanism through 
two energetically closely spaced molecular states. The first 
state (red) couples only to the majority spin band, while the 
minority spin couple only to the purple state. In the parallel 
case one finds T n (E) / and T U (E) / but in the antipar- 
allel T n (E) = 0. We then expect an infinite GMR for such 
an energy window. 



demonstrating that the GMR signal can be tuned by an 
appropriate choice of the anchoring chemistry. 
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Figure 15: Transport properties for a 1,4-phenyl molecule 
attached to Ni (100) surfaces through a S group. The top 
panel shows the I-V characteristics for both the parallel and 
antiparallel alignment of the leads and the inset the corre- 
sponding GMR ratio. The lower panel is the transmission 
coefficient at zero bias as a function of energy. Because of 
spin-symmetry, in the antiparallel case we plot only the ma- 
jority spin. Reprinted with permission from— A.R. Rocha et 
al., Phys. Rev. B 73, 085414 (2006). Copyright American 
Physical Society 2006. 



Finally it is worth mentioning that a severe bias de- 
pendance of the GMR ratio, with also the possibility of 
negative values, has been recently predicted with either 
semi-empirical^ and ab initial methods. 
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Figure 16: Transport properties for a 1,4-phenyl molecule at- 
tached to Ni (100) surfaces through a Se group. The top 
panel shows the I-V characteristics for both the parallel and 
antiparallel alignment of the leads and the inset the corre- 
sponding GMR ratio. The lower panel is the transmission 
coefficient at zero bias as a function of energy. Because of 
spin-symmetry, in the antiparallel case we plot only the ma- 
jority spin. Reprinted with permission from— < A.R. Rocha et 
al, Phys. Rev. B 73, 085414 (2006). Copyright American 
Physical Society 2006. 



B. Carbon Nanotubes 

Carbon nanotubes are almost defect-free graphene 
sheets rolled up to form one-dimensional molecules with 
enormous aspect ratios*^. Their conducting state (metal- 
icity) depends on their chirality, however in the metallic 
configuration they are ideal conductors with a remark- 
ably long phase-coherence length—^. An important as- 
pect is that the relevant physics at the Fermi level is 
entirely dominated by the p z orbitals, which are radi- 
ally aligned with respect to the tube axis. These in- 
clude the bonding properties with other materials and 
between tubes. Therefore carbon nanotubes appear as an 
ideal playground for investigating both GMR and TMR 
through molecules. In fact one can expect that two tubes 
with different chirality will bond to a magnetic surface in 
a similar way, allowing us to isolate the effects of the 
molecule from that of the contacts. Indeed TMR-like 
transport through carbon nanotubes has been experi- 
mentally reported by several groupo 20 i 79 ' 80 i 81 ' S2 i 83 ' 84 i 85 . 

Why would one expect a large GMR from a carbon 
nanotube? To answer this question we use an argument 
derived by TersoflS^ and then subsequently refined 8 ^ for 
explaining the contact resistance between a C-nanotube 
and an ordinary metal. Consider for simplicity an arm- 
chair nanotube (metallic). The Fermi surface of such 
tube consists only of two points, symmetric with respect 
to r in the ID Brillouin zone (see figure IT8|) . The Fermi 
wave-vector is then ftp — 2tt/3zq with zq — o?oV / 3/2 and 
d the C-C bond distance (d —lA2 A). 



Figure 17: Transport properties for a 1,4-phenyl molecule at- 
tached to Ni (100) surfaces through a Te group. The top 
panel shows the I-V characteristics for both the parallel and 
antiparallel alignment of the leads and the inset the corre- 
sponding GMR ratio. The lower panel is the transmission 
coefficient at zero bias as a function of energy. Because of 
spin-symmetry, in the antiparallel case we plot only the ma- 
jority spin. Reprinted with permission from— A.R. Rocha et 
al, Phys. Rev. B 73, 085414 (2006). Copyright American 
Physical Society 2006. 
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Figure 18: Fermi surfaces of an armchair carbon nanotube 
and of a magnetic transition metal. The Fermi surface of 
the nanotube consists in two points fcp = q, symmetric with 
respect to the Y point. The Fermi surface of a transition 
magnetic metal consists of two spheres (for j and [ spins) 
whose different diameters depend on the exchange field. The 
three possible scenarios discussed in the text: (a) q < fcp < 
4, (b) 4 < 4 < q, (c) 4 < 9 < 4- 

Assume now for simplicity that our magnetic metal is 
an exchange-split free-electron gas, whose band-energy is 

h 2 k 2 

E° k = — + aA/2, (12) 
2m 

with a = — 1 (a = +1) for majority (minority) spins and 
A the exchange energy. The spin-dependent Fermi wave- 
vectors are then respectively k F = ^ s /2m(E-p + A/2)/fr 

and 4 = V 2m ( E F - A/2)/n. 

The transport through an interface between such a 
magnetic metal and the nanotube is determined by the 
overlap between the corresponding Fermi surfaces. Three 
scenarios are possible. First the Fermi- wave vector of the 
carbon nanotube is smaller than both 4 an d 4 ( see -^S" 
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ure 118b ) . In this case in the magnetic metal there is 
always a ft-vector that matches the Fermi-wave vector of 
the nanotube for both spins. Therefore both spins can 
be injected into the tube and the total resistance will be 
small and spin- independent. 

Secondly the Fermi-wave vector of the carbon nan- 
otube is larger than both fcX and ftp (see figure 118b ) . 
Now no states are available in the metallic contacts whose 
wave- vectors match the Fermi wave- vector of the carbon 
nanotube. For zero-bias and zero-temperature the re- 
sistance is extremely large. Nevertheless upon increasing 
the temperature, phonon assisted transport becomes pos- 
sible. Spin electrons can be scattered out of the Fermi 
surface into states with large longitudinal momentum. 
At temperature T the fraction of electrons with energy 
above E-p is simply proportional to the Fermi distribution 
function. However, because of the exchange energy, spin- 
up electrons will possess higher momentum than spin- 
down. Therefore one can find more spin-up states with 
a longitudinal momentum matching the one of the nan- 
otube than spin-down states. This potentially gives a 
temperature-induced spin-dependent resistance. 

Finally if the Fermi wave-vector of the carbon nan- 
otube is larger than ftp but smaller than fc F (see figure 
118b ) . only the majority electrons can enter the nanotube. 
The system becomes fully spin-polarized and the Fermi 
surface matching replaces the bonding spin-selectivity en- 
countered in the previous section. Also in this case a spin- 
valve structure made by magnetic contacts and carbon 
nanotube as spacer is predicted to show infinite GMR at 
zero temperature. An increase of the temperature will 
produce a degradation of the polarization since minor- 
ity spins can be thermally scattered outside their Fermi 
surface, and therefore contribute to the transport. 

Two important aspects must be pointed out. First all 
these considerations are based on the assumption of per- 
fectly crystalline systems. This may not be true in reality 
and the effects of breaking the translational invariance 
must be considered. From a qualitative point of view 
disorder smears the Fermi surface and eventually may 
produce some states with large longitudinal momentum. 
This will improve the conductance through the nanotube, 
even if its spin-polarization will be in general dependent 
on the nature of disorder. 

Secondly, our heuristic argument does not necessar- 
ily apply to the case of transition metal contacts, where 
the spin-selectivity arising from the different bonding 
nature of the two spin-subands can play an important 
role. Clearly more realistic bandstructure calculations 
are needed. These however are rather problematic. In 
addition to the need of simulating transition metal leads 
one has to consider rather large supercell for a reasonable 
description of the nanotube/ metal interface. 

For this reason most of the calculations to date have 
used simple tight-binding models without self-consistent 
proceduresSi22i2S. These roughly agree on the possibility 
of large GMR ratios in transition metals contacted nan- 
otubes, although the actual values predicted are some- 



how affected by the different methods and the contact 
geometry. 



C. Long molecules and phonons 

As mentioned previously electronic transport in poly- 
mers and long molecules offer a considerable higher 
level of complexity when compared to small ballistic 
molecules. The first consideration is that the trans- 
port is often driven by strong electron-phonon interac- 
tion, therefore the relevant current carriers are some cor- 
related electronic-vibronic states (polarons, bi-polarons, 
etc.). Moreover in the case of spin-transport one has 
to take into account the likely presence of paramagnetic 
centers^. They usually appear in low concentrations, 
but they become relevant to the spin-dynamics in long 
molecules. These two features add to the large scale of 
the system and make the problem of spin-transport in 
large organic molecules intractable with ab initio meth- 
ods. This is the main reason why to date only simple 
Hamiltonian models have been used. 

When the electron-phonon interaction is weak, or al- 
ternatively the molecule is rather short in such a way that 
polaronic-type of transport is not dominant, then phonon 
absorption/emission has only the effect of smearing the 
transmission coefficient. In spin- valves made from transi- 
tion metals such a smearing is likely to result in a reduc- 
tion of the spin-polarization of the device and therefore 
of the GMR 91 . However quantitative predictions are dif- 
ficult, since also in this case a detailed description of the 
electronic structure of the electrodes and of the bond- 
ing with the molecule is essential. For instance calcula- 
tions of spin-injection into short strands of DNA obtained 
with simple single-orbital tight-binding models and some 
parameterization of the metal/molecule bonding, report 
both an enhancement^ and a reduction 9 ^ of the GMR 
with the bias. 

When the electron-phonon interaction is strong and 
the molecules are rather long, then the ground state of 
the molecule is some correlated electronic-vibronic state. 
In this case the situation is more complex and to the best 
of our knowledge no transport calculations have been car- 
ried out to date. Some interesting insights come from the 
investigation of the ground state of polymers in contact 
with magnetic metals. Xie and co- workers^ investigated 
a non-degenerate polymer in contact with a model metal 
reproducing either a magnetic transition metal or man- 
ganite. Interestingly they found that when no charge is 
transferred from the metal to the polymer, this remains 
metallic with a rather uniform distribution of the charge 
across the interface. In contrast charge transfer promotes 
the formation of spinless bipolarons in the polymers. An 
analysis of the DOS of the whole structure further sug- 
gest that the bipolaron formation drastically reduces the 
spin polarization of the whole system. 
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VI. MORE EXOTIC PHENOMENA 
A. Molecular Magnets 

Extremely interesting features appear in quantum 
transport through molecules when the conducting elec- 
trons interact with some internal molecular degrees of 
freedom. This is for instance the case of vibrational 
levels, where electron-phonon interaction manifests itself 
with sharp changes in the slope of the I-V curve, steps in 
dl/dV and in peaks in d 2 I/dV 2 . It is therefore natural 
to speculate about similar effects associated to internal 
magnetic degrees of freedom, i.e. the molecular spin 94 . 
Until recently the experimental activity was focussed on 
reproducing at the molecular level effects already demon- 
strated in quantum dots, such as the Zeeman splitting of 
the Coulomb-blockade as well as the Kondo effec t 95 ' 96 . 
More recently, advances in the chemical functionaliza- 
tion of magnetic molecules have allowed the formation 
of stable bonding between the molecules and metallic 
surfaces2L2£, thus the construction of two and three ter- 
minal devices22ii2S. 

Magnetic moleculesifli are molecules comprising a 
number of transition metal ions magnetically coupled to 
each other in such a way to give rise to a global net 
spin S. They are usually described by the following spin 
Hamiltonian Hq 

H =DS 2 +g^ B H z S z , (13) 

where D (D < 0) is the zero field splitting constant and 
H z is the z component of an external magnetic field (/ib 
is the Bohr magneton). The term proportional to S 2 
lift the degeneracy of the spin-multiplet and the energy 
levels can be labeled by the magnetic quantum number 
M s , with -S < M s < S. 

If the dynamics is solely determined by the Hamilto- 
nian Hq then a molecule prepared in a given Ms state 
will remain in such state. However a perturbation H±, 
which does not commute with Hq will create mixing be- 
tween the Ms states, thus transition between levels with 
different Ms will be possible. Selection rules for these 
transitions are given by the particular symmetry of the 
perturbation, and for instance transverse anisotropy 

ffi = E{S 2 X - S 2 ) , (14) 

allows transitions Ms — > Ms ± 2n with n an integer. 
Evidence of quantum tunneling of the magnetization for 
magnetic molecules are now numerous^fli. 

Let us now consider the case of electronic transport 
through such a molecule and in particular the case 
of weak coupling between the molecule and the cur- 
rent/voltage electrodes, which is the most likely situa- 
tion in actual experiments22ii2S. In general the ground 
state for a charged molecule will be different from that 
of its neutral state, and so will be the excitation spec- 
trum. This means that sequential tunneling can poten- 
tially leave the molecule in some excited state not allowed 



to relax by the selection rules. Importantly also the op- 
posite is true, namely that a charged exited molecule can 
relax to a state which does not allow electron transfer to 
the electrodes because of the selection rules. This may 
lead to a complete suppression of the currents^ 9 ". 

Two experimental works have been published so 
far on transport through magnetic molecules 99,10 -. 
They both consist in a transistor geometry 
where a single Mni2 molecule is trapped. Mn.12, 
[Mni 2 Oi2(CH 3 COO)i 6 (H20)4], is perhaps the proto- 
type of all molecular magnets. The 12 Mn ions occupy 
three inequivalent atomic sites, namely two Mn 3+ and 
one Mn 4+ . The Mn ions with different valence couple 
antifcrromagnetically to each other resulting in a total 
S = 10 ground state. In one case the molecule has been 
functionalized (see figure I19f) with thiol groups in order 
to form stable contacts with the gold surfaces of the 
electrodes. The crucial aspect of both these experiments 
is the discovery of fingerprints of the molecular state 
into the I-V characteristics of the device. In particular 
negative differential conductance, complete current 
suppression^ and non-trivial dependence of the peaks 
in the differential conductance over an applied magnetic 
fieldiffi, can all be interpreted as a result of the internal 
spin-dynamics of the molecule. 




Figure 19: Top and side views of a ball-and-stick model 
for a Mni2 magnetic molecule. In this case we con- 
sider [Mni20i2(0 2 C-R-SAc)i6(H 2 0)4] where R={ C 6 H 4 , 
C15H30 }. Colour code: green=carbon, yellow=sulphur, 
blue=hydrogen, red=oxygen and violet = manganese. 

Several model calculations have been performed for 
explaining the various features of these experiments^^ 
and proposing a new setup where one of the electrodes 
is a ferromagne1ii22*i2i. These are based on the spin- 
Hamiltonian of equations 113|) and l|14|) and an addi- 
tional term that takes into account the coupling with the 
electrodes. Then the transport calculation is performed 
within the standard master equation formalism familiar 
to transport through interacting quantum dotsiS 5 -. Al- 
though these calculations are certainly important since 
they demonstrate the mechanisms behind the various ef- 
fects, they depend heavily on the specific choice of the 
parameters used. For this reason ah initio simulations 
would be very important. For instance questions about 
the actual charging state of the molecule under bias, the 
strength of the interaction between the conducting elec- 
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trons and the Mn ions, and the electrical response of the 
whole system are likely to find an answer with DFT cal- 
culations. The problem is indeed complicated and even 
the simple evaluation of the ground state of magnetic 
molecules from DFT is not triviali2Sii2i. For this reason 
no first principle calculations of transport through mag- 
netic molecules have been performed to date. However we 
believe that this is an extremely challenging field where 
order-N capability, strong correlation and scalable quan- 
tum transport schemes can find a common playground. 

Finally we want to discuss the possibility of construct- 
ing all-molecular spin- valves, i.e. spin- valves where both 
the spin-injector and spin-detector are molecules or part 
of the same molecule and the metallic current /voltage 
electrodes have the only function of electron reservoirs. 
Also in this case the expectation is to detect electrically 
the magnetic state of the molecule, however for spin- valve 
operation we also demand the possibility of switching re- 
versibly between two spin-configurations by applying a 
magnetic field. An interesting proposal is that of using 
dicobaltocene moleculesi2&. These belong to the metal- 
locene family and are characterized by two Co ions sep- 
arated by a spacer that can be chemically engineered. 

DFT calculations^ show that dicobaltocene attached 
to gold electrodes in its ground state is stabilized by su- 
perexchange in an antiferromagnetic configuration, i.e. 
the local moments of the two Co ions are antiparallel to 
each other. As in a conventional spin- valve a parallel 
alignment can be obtained by applying a large magnetic 
field (around 20 T for a C2H4 spacer), and the calcula- 
tion shows a rather large GMR ratio. This is suggestive 
of the possibility of all-molecular spin- valves and future 
experiments in this area are certainly welcome. 



B. d ferromagnetism and magnetic proximity 

We wish to close this review with a brief discussion of 
two recently discovered intriguing phenomena, that chal- 
lenge our current understanding of ferromagnetism and 
may offer a new playground for spin-transport. These are 
d° ferromagnetism and magnetic proximity effect. The 
measurement of long-range ferromagnetic order in ma- 
terials not containing ions with either d or / electrons, 
therefore without an obvious way to produce a local mag- 
netic moment, is the common factor of these two aspects 
of magnetism. 

Let us consider first the case of d° ferromagnetismiSS. 
This is the intrinsic ferromagnetism of usually highly de- 
fective materials, which do not contain any partially filled 
d or / shells. Among the several examples we wish to 
mention irradiated graphite and fullerenesiiS, nonstoi- 
chiometric CaBgiii, and Hf02 thin filmgAi. 2 . A com- 
mon explanation for the magnetism in all these mate- 
rials is not available at present. Note that one has to 
explain both the formation of a magnetic moment and 
the long-range coupling between moments. Generally the 
moments are associated to intrinsic defects. Paramag- 



netic defects in organic materials are not uncommon 60 
and strongly correlated molecular orbitals associated to 
vacancies have been suggested for magnetic oxide a 113 i 114 . 
However the demonstration of long range coupling re- 
mains elusive. For instance recent DFT calculations^^ 
demonstrate, at least for CaO, that the defect concentra- 
tion needed for long range ferromagnetism is three orders 
of magnitude larger than that obtainable at equilibrium. 
Finally we wish to mention that ferromagnetism origi- 
nating from p-shells has been suggested for oxygen mixed 
valence compounds such as Rb40gii&. 

In contrast to d° ferromagnetism magnetic proximity 
is not an intrinsic material property and it does require 
the presence of a magnetic material. The basic idea is 
quite simple: there is always some charge transfer at the 
contact between a conducting molecule and a metal as- 
sociated with the alignment of their respective chemical 
potentials. In a ferromagnet some degree of spin trans- 
fer accompanies the charge transfer giving rise to an in- 
duced magnetic moment. Note that the magnetic prox- 
imity effect does not imply intrinsic ferromagnetism and 
no spin aligning potential exists in the non-magnetic ma- 
terial. This means that a magnetic moment is detected 
only when a second ferromagnetic material is present and 
when good contact is made. 

A first indirect evidence of this effect was found in 
explaining the unaccounted magnetization in a carbon 
meteorite rich of magnetic inclusionsiit Then Ferreira 
and Sanvito derived a close system of equations for the 
induced magnetic moment and for the energetic of a car- 
bon nanotube deposited on a magnetic surface^. The 
calculation, based on a simple tight-binding model, re- 
vealed that induced magnetic moments of the order of 
0.1 /ib per carbon atom in contact can be achieved at 
room temperature. This paved the way for a more con- 
trolled set of experiments. 

Experimentally the problem is to detect the small in- 
duced moment over a huge background coming from the 
ferromagnetic substrate. One possible strategy is that of 
measuring the stray field around the nanotube once this 
is placed on a smooth thin film. A uniformly magnetized 
thin film creates no stray field whatever its direction of 
magnetization. In contrast a magnetized nanotube will 
produce a stray field, which will be directly detectable. 
This idea was exploited in the experiments from Cespedes 
et al^iS, who measured induced magnetic moments in ex- 
cess of 0.1/zb per carbon atom in contact, in good agree- 
ment with the theoretical prediction. The experiment 
consists in taking AFM and MFM images of nanotubes 
on various surfaces. The difference between the topo- 
graphic (AFM) and the magnetic images (MFM) is a di- 
rect measure of the stray field coming from the nanotube 
and therefore provides evidence for the induced magnetic 
moment (see figures 12171 and |2*T1) . 

In the another experiment, Mertins et alm^ produced 
a multilayer of thin, alternating iron and carbon layers, 
of thickness 2.55 and 0.55 nm, respectively. Then, they 
probed locally the magnetic moment of the carbon by X- 
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Figure 20: AFM images of a carbon nanotube on copper a) 
and silicon c) substrate and their corresponding MFM scans 
b) and d). The MFM images do not show any magnetic con- 
trast indicating no induced magnetic moment. 

ray magneto-optical reflectivity of polarized synchrotron 
radiation. In this type of measurement the Fe and C 
absorption edges differ by about 500 eV enabling one 
of establishing with precision whether the magnetic mo- 
ment comes from C or not. With this method magnetic 
moments of the order of 0.05 /zb were found. 



Sllil 




Figure 21: AFM images of a carbon nanotube on cobalt sub- 
strate (left) and the corresponding MFM scans (right). The 
MFM image shows magnetic contrast indicating an induced 
magnetic moment. 

Finally we wish to mention a few experiments where 
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VII. CONCLUSIONS 

We have reviewed the most recent advances in spin- 
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ticated algorithms capable of high accuracy and of un- 
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